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We study the interplay between rotating wave approximation and optimal control. In particular, 
we show that for a wide class of optimal control problems one can choose the control field such that 
the Hamiltonian becomes time-independent under the rotating wave approximation. Thus, we show 
how to recast the functional minimization defined by the optimal control problem into a simpler 
multi-variable function minimization. We provide the analytic solution to the state-to-state transfer 
of the paradigmatic two-level system and to the more general star configuration of an A'^-level system. 
We demonstrate numerically the usefulness of this approach in the more general class of connected 
acyclic A'^-level systems with random spectra. Finally, we use it to design a protocol to entangle 
Rydberg via constant laser pulses atoms in an experimentally relevant range of parameters. 


The Rotating Wave Approximation (RWA) plays a ma¬ 
jor role in simplifying the quantum mechanical descrip¬ 
tion of laser driven systems in the low intensity regime: 
it takes into account only the co-rotating field with the 
system and it neglects the counter-rotating part fli- 
This approximation has been introduced for two level 
quantum systems, and then generalized for A^-level sys¬ 
tems 0 . The deviations from this approximation for big 
intensities are well known and commonly described as 
Bloch-Siegert shifts, breaking the harmonicity of the sys¬ 
tem dynamics Q. Finally, a more general description is 
given by Floquet theory that allows to treat periodically 
driven systems 

Developing error-free protocols for the manipulation of 
quantum systems - also along the development of quan¬ 
tum technologies but not restricted to them - is one of 
the major challenges in contemporary research in atom 
and molecular physics During the last decades, an 
increasing contribution in such effort has come from the 
exploitation of Quantum Optimal Control (QOC), the 
search for an optimal control pulse to perform a given sys¬ 
tem manipulation Q. Methods to solve QOC problems 
have been developed [i®] and experiments have shown 
the great benefit from them, see e.g. [iMli- We have 
now deep theoretical understanding of QOC, in particu¬ 
lar about the Mssibilities and hurdles to control quan¬ 
tum systems |16l-lla|. and we even start to understand 
the complexity of QOC problems [H, Graph the¬ 
ory concepts have been exploited to attack a question 
that lies at the heart of controllability studies: given a 
Hamiltonian depending on some time dependent tunable 
control field, is it possible to dynamically connect every 
pair of arbitrary initial and final state, i.e. it is possi¬ 
ble to realize every possible state-to-state transfer? A 
widely used criterion to answer to this question is via 
dynamical Lie-Algebras [l^, while here we make use of 
a different criterion based on graph theory: Turinici and 
Rabitz showed that if the graph corresponding to the 
control Hamiltonian is connected and the spacings of the 
eigenvalues of the uncontrolled part of the Hamiltonian 
are non-degenerate, then wave function controllability is 
guaranteed El- 


In this paper, we identify a significant class of QOC 
problems where the RWA is applicable and show that 
within this class QOC problems can be easily solved: 
the functional optimization can be recast into a multi- 
variable function minimization, thus simplifying the nu¬ 
merical efforts and improving the theoretical understand¬ 
ing of the process. For this purpose we employ concepts 
from graph theory to analyze the system dynamics in an 
easier and more practical picture [2lj. This allows us 
to straightforwardly identify the cases where it is possi¬ 
ble to map the time-dependent Hamiltonian onto a time- 
independent one and to find the solution of the QOC 
problem. As a result of this analysis, we can show that 
in this class of QOC problems the lower bound for the 
number of independent control parameters necessary for 
a successful control of the system - introduced in a recent 
work - is saturated [l^ . 


Finally, as an example of possible applications of this 
approach we use it to design a protocol to entangle 
Rydberg atoms [l^. Indeed, in the last years experi¬ 
ments with Rydberg atoms have attracted increasing at¬ 
tention as a promising platform for implementing quan¬ 
tum information processing alg orithms, such as CNOT- 
gates [IJ , two-body [2§-[27 1 or many-body entangle¬ 

ment 12811. 


The paper is structured as follow: In Section U we 
present the model studied and for the sake of complete¬ 
ness we review the work of Einwohner et al. on the gen¬ 
eralized Al-level RWA Q and in Section |H] we review the 
algorithm based on graph theory to recast systems in a 
time-independent form. Section fllll presents the analytic 
solutions of the state-to-state transfer in the two-level 
system and a special case of the Al-level system while 
numerical results are presented in Section El In Sec¬ 
tion El we demonstrate the methods presented here by 
designing a protocol entangling Rydberg atoms. Finally, 
in Section IVll we discuss the results and present the con¬ 
clusions of this work. 
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I. THE MODEL 

In the standard RWA setting one considers a system 
described by a A^-dimensional Hamiltonian {h = 1) 

F 

/=i 

where Ho is the time-independent drift Hamiltonian 
with eigenvectors \n) and corresponding eigenvalues 
En in = 0,1,...A^ — 1) and He is the control Hamil¬ 
tonian, which might correspond to different physical sce¬ 
narios. Here, we only assume that the diagonal elements 
of He are all zero. The coefficients Af G C represent the 
independent control parameters and w/ G R the driv¬ 
ing frequencies. A relevant example is naturally encoun¬ 
tered in molecular or atomic physics when describing the 
interaction between matter and light in the dipole ap¬ 
proximation, where each control with frequency ujf and 
strength Af is typically realized by a laser and He is 
the dipole operator. In the following we assume that 
no two driving frequencies are in approximate resonance 
to the same transition frequency E^j = Ek — Ej and on 
the other hand every driving frequency is in approximate 
resonance to at least one transition frequency. There is 
not much loss of generality in this when considering a 
QOC problem since this only means that we do not have 
off resonant controls nor two controls affecting the same 
transition. Assuming that the transition frequencies are 
not degenerate {Ekj ^ Ek'f V (fc,j) A (k\ j')) we can 
make use of results presented in Ref. [l7| that ensure 
wave function controllability. The last assumption we 
make is that the RWA is valid, i.e. we are in a setting 
of low intensities and resonant driving frequencies. That 
is we assume for each set of transition frequency Ekj, 
with driving frequency w/ in approximate resonance to 
that transition, and the corresponding amplitude Af the 
inequalities 

LOf ^ Akj = \Ekj — ujf\ 

to hold. 

Expanding the wave function \tpit)) as 

N-l 

\m) = Y,Ckit)e-^^'‘*\k), (3) 

the Schrodinger equation 

^±\i;{t))=Hm)) ( 4 ) 

can be rewritten as a differential equation for the coeffi¬ 
cients c = (cq. Cl,..., cvr_i): 

N-l 

ick{t) = (5) 

1=0 


where 

F 

= 2 ^ + A}Ei^’‘^+‘^fi\He)ko (6) 

/=i 

Following Q , the multilevel rotating wave approximation 
is done by neglecting all non-resonant terms H^^. This 
includes far more than just dropping the counter-rotating 
term (one of the terms inside the sum of Eq. (| 6 |)), as in 
the more common RWA for two-level systems, but also 
all terms of the sum of frequencies that are non-resonant 
to the transition frequency corresponding to that matrix 
element, reducing the whole sum to only one element. 
The resulting matrix elements are denoted by Hj^p = 
, with 

^ [(1 + sign(Efej))^/ + (1 - sign{Ekj))A}] . 

( 7 ) 

A change of basis bk{t) = e^^’‘*Ck{t), G R which has no 
a priori physical meaning but is merely a mathematical 
tool transforms Eq. m into 

N-l 

ibk{t) = Y ( 8 ) 

i=o 

To end up with a time-independent operator one has to 
set all phases equal to zero, that is all ^(A^— 1 ) equations 

ilk - H =0, (fc ^ j) (9) 

have to be solved by choosing an appropriate 7 = 
( 71 , 72 ,..., 7 fc). The resulting system has the solution; 

hit) = exp i — diag( 7 )^ = 0 ). ( 10 ) 

The term iHc)kj is included in ([9]) since only for non¬ 
vanishing matrix elements the corresponding equation 
has to be considered. In general those equations ([9]) can¬ 
not be solved consistently, but the zeros in (i?c)fcj de¬ 
crease the number of equations to be solved. This trans¬ 
formation has the advantage of possibly making the de¬ 
scription of the system simpler. In particular for acyclic 
graphs we will see that system can always be reduced to 
a time-independent one. 


II. TIME-INDEPENDENT DESCRIPTION FOR 
ACYCLIC GRAPHS 

Every Hamiltonian consisting of a diagonal drift term 
Ho and a control term He has a pictorial representation 
via its energy level scheme, which can easily be mapped 
onto an undirected graph, see Fig. [T] for a descriptive 
example. Every state {Ej} is portrayed by a vertex and 
edges of the graph illustrate non-zero transition elements. 
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them. Repeating this step, that is removing i?i, leaves 
us with the root graph: Eq- In the last three pictures 
we show how the variables 7fc are chosen. This is done 
by assigning an arbitrary real weight 70 to Eq and then 
re-adding the pendant vertices Ei , E2 and E4 with their 
corresponding weights 7fc = 7 — Afco,fc = 1 , 2 , 4 . The 
algorithm finishes with adding the last vertex E3 with 
its weight 73 = 71 — A31. 


FIG. 1: (Color online) Energy level scheme of an example 
five-level system (left-hand side) and the corresponding graph 
(right-hand side). The graph is connected but has a cycle 
(Eo, Ei,Eo). Removing the red dotted edge between the ver¬ 
tices Eo and Eo makes the graph acyclic while maintaining 
the connected property. Removing the blue dashed edge be¬ 
tween El and Eo also leaves the graph connected and acyclic, 
but in addition that graph is a star, a type of graph which is 
studied in Sec. niBi 


There is a specific class of graphs (and therefore control 
Hamiltonians), that deserve a special treatment due to 
their convenient properties. Those are the subset of con¬ 
nected, acyclic graphs (“trees”), where from every vertex 
k to every other vertex j there exists a unique path, so 
no cycles are present. In Fig. [T] removing the blue dashed 
edge between Ei and E3 or the red dotted edge between 
Eq and Eq from the complete graph results in an acyclic 
graph. As already mentioned in the introduction, con¬ 
nected, acyclic graphs form a very significant subset of 
QOC Hamiltonians: They connect the complete system 
dynamics with the smallest number of control parame¬ 
ters: all complex amplitudes Af and frequencies w/ are 
independent parameters to control. As a consequence the 
number of controls is set to F’ = — 1 and the number of 

equations to solve reduces to A — 1, so a solution for the 
equations m always exist. In the following we briefly ex¬ 
plain a constructive approach to find a solution, referring 
again to Ref. Q : 

1 . Every pendant vertex has a unique “successor” ver¬ 
tex, connected to it by an edge 

2 . Delete recursively all pendant vertices until one 
ends up with one vertex I, assign an arbitrary real 
value 7; to it. 

3 . Every vertex k is assigned a value jk = Ij — ^kj, 
where j is the successor of k. 

The algorithm will be implicitly used in Sec. IHI Al and 
IIII Bl and is directly implemented in the numerical calcu¬ 
lations presented in Sec. IHI HI 

To elucidate this algorithm we depicted an example 
based on the graph from Fig. [T] without the red dotted 
edge (see Fig. [ 2 |). The first three pictures illustrate the 
“pruning” of the tree, that is finding the root of the tree. 
That is done by finding all pendant vertices of the graph, 
i.e. those that are only connected via one edge, which 
are in this case the vertices E2 , Eq and Ei , and removing 


III. ANALYTICAL RESULTS 

In the following we present two systems that have a 
simple analytic solution: the paradigmatic two-level sys¬ 
tem and a A-level system whose graph is a star. Those 
solutions provide a general insight in the structure of dy¬ 
namics and hence the possibilities for QOC. The algo¬ 
rithm reviewed in Sec. m is implicitly used in the follow¬ 
ing parts. 


A. State-to-state transfer for the two-level system 

As a starting point of our analysis, we specialize our 
investigation to the case where A = 2 and consequently 
F = 1 . In this case the system given by © takes the 
form 


H = HD+Re{AE‘^*)Hc, ( 11 ) 

Ac = |0)(l| + |l)(0|, (12) 

where the indices of A = Ai, w = wi and A = Aqi were 
dropped for clarity. Employing the RWA and assuming 
the system initially to be in the ground state c (t = 0) = 
(1,0), one recovers the known solution 

A A A \ 

cos —t — i— sm — t 
2 A 2 J 

ci(t) = e-*^/^‘ 4 sin-t, 

^ ^ A 2 ’ 

where A = ^A'^ + |Ap. Using the Bloch sphere 
representation and cg( 0,(/)) = (cos sin ^ G 

[O, 27 r ],0 G [0,7r] as the goal state for a state-to-state 
transfer, we find the equations 



( 13 ) 

( 14 ) 


,A/2/ A A . A \ 0 

2 ' cos—t —z—sm—t =cos — 

I 2 A 2 2 


( 15 ) 


e-*^/2*^sin-t = sin-e*'^. ( 16 ) 

A 2 2 ^ ^ 

From Eq. (HU) one can find that for a given goal state 
the tunable parameters of the system have to satisfy the 
inequality 


2 arcsin -LJ. > 0 
A “ 


( 17 ) 
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FIG. 2 : Graph corresponding to an energy scheme of a five-level system taken from Fig. [T] except for the red dotted edge. The 
first three figures show the “pruning” of the tree, where all pendant vertices and the edges connecting to them are deleted 
recursively until one root (So) is left. The last three figures depict how the weights {7*:} are recursively assigned to the edges: 
First the root So is assigned an arbitrary value 70, then the three neighboring vertices Si,S2,S4 are assigned their values 
7J, = 70 — Afco, A: = 1 , 2 , 4 and at last one assigns S3 its value 73 = 71 — A31. 


and that the total time T is 

T = 2 ^rcsin(Asin|). (18) 

This time T is known in the literature as the quantum 
speed limit, that is the smallest time necessary to evolve 
in Hilbert space from the initial state to the goal state at 
a given fixed energy [ 2 ^. Using Eq. m and comparing 
the complex phases we see that the phase a of A = |H|e®“ 
has to be chosen as 


a = (j)+ (19) 

In summary the three Eqs. (IT71) . (fTBl) and (fTOl) yield all 
information necessary to control the system: Given a pair 
(|H|, A), (I17|) defines the Bloch vector with the maximal 
distance on the Bloch sphere from the initial state that 
can still be reached, thus every state with an angle 0 
smaller than this value can also be reached. Those states 
define the set of reachable angles. Note that for A = 0 
we can reach all states on the Bloch sphere. The time 
specified in Eq. (ITO)) provides the time necessary to per¬ 
form such a process if the desired Bloch vector is exactly 
at the boundary of the set of reachable states given by 
Eq. (IT71) . Finally, Eqs. (IT51) and (fTOl) together provide 
conditions for the amplitude and phase of the control. 


B. State-to-state transfer for the star A-level 
system 

For a generic TV-level system, even with an acyclic 
graph, analytic solutions are rare. However, if we restrict 
the interaction furthermore such that only the transitions 
| 0 ) |fc) (fc 7 ^ 0 ) are allowed and increasing the symme¬ 
try in the problem by setting all detunings to an equal 
value (Aofc = A Vfc 0), an analytic solution can be 
found. In the example presented in Fig. [T] this scenario 
occurs if from the complete graph the blue dashed edge 
between Ei and E 4 is erased. 

For a clearer display of the underlying dynamics we again 


assume the system initially to be in the ground state 
c{t = 0) = (1,0,0,..., 0) of the drift Hamiltonian and 
define analogously to Sec. IHI Al 




\ 


N-1 

/=i 


( 20 ) 


After some straightforward algebra, one finds 
co(t) = “ 

Cfc(t) = e“*^/^*^sin^t (fc = 1, 2,..., TV - 1). (22) 
z 


A A \ 


To solve the state-to-state transfer problem 
we introduce the normalized goal state ^ = 

where and Pk 
are real numbers and the global phase is chosen such 
that Po = 0. The procedure from here on and the 
structure of the resulting equations follow that of the 
two-level system: for a goal state to be reachable, the 
inequalities 


^ > Cfc (23) 

have to hold for all fc 7 ^ 0. The total time T is then fixed 
by 

T = max arcsin ^ ^ ^ 


Here, T is again what is referred in the literature as the 
quantum speed limit. Furthermore the phases ak of Ak = 
|Afc|e®“'“ have to be chosen as 

ak=Pk + ^AT. (25) 

As before, the Eqs. (1^ . (IMl) and (E5|) determine if and 
how a goal state can be reached. While Eq. (1^ refers 
to the reachability, Eqs. (1^ and (E5|) state how the am¬ 
plitudes and phases of the different Ak have to be chosen 
to solve the problem provided that the goal state is at 
the boundary of the reachable set. Note that again for 
A = 0 we can reach all states. 
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IV. NUMERICAL RESULTS 

To bolster the analytic results obtained so far and 
show the frontier up to which the RWA is an excel¬ 
lent tool for QOC we simulate systems discusses in 
Sec. El We create random non-degenerate energy spec¬ 
tra {Aq, Ai,..., A^v-i}, connected and acyclic Control 
Hamiltonians He and use the algorithm mentioned in 
Sec. El to set up a time-independent Hamiltonian. As a 
figure of merit for state-to-state transfer we use the in¬ 
fidelity J = 1 — |('0(3|'0(r))|^. The parameters for the 
optimization are the amplitudes |Ay|, the phases af and 
the final time T. We fix the detuning A (which is as¬ 
sumed to be equal for all frequencies to study the effects 
of the detuning via changing one parameter instead of an 
increasing number while adding more freedom the task 
of optimization would be easier) and optimize for many 
different random goal states. The optimization is per¬ 
formed via the direct search method Nelder-Mead (3nl |. 
To check if the results obtained in the RWA are valid in 
the complete description or if the RWA breaks down, we 
perform an exact time evolution with the results of the 
optimization and compared the results. Again, the infi¬ 
delity is used as a figure of merit. 



FIG. 3: (Color online) Percentage A of successfully optimized 
states as a function of detuning A for different dimension 
N: 2 (green), 3 (blue), 4 (pink), 5 (cyan) and 6 (yellow). 
Inset: Percentage A of positive double check as a function 
of detuning A for N = 2, 3, 4, 5 and 6 with the same color 
scheme as the main figure. 


In Fig. [3] the percentage of randomly generated states 
reached A within an infidelity of less than e = 10“^ 
is shown as function of the detuning. We see that it 
increases with decreasing detuning and that for every 
detuning smaller than 10“"^ we are able to reach all 
states, that is the percentage of reached states is equal to 
A = 100%. Moreover, for bigger dimensions the percent¬ 
age of reached states decreases if the detuning is bigger 
than 10“"^. The inset shows the result of the numeric 
double check, where we computed the time evolution of 


the exact system (without the RWA) and then calcu¬ 
lated the infidelity of the exact time evolved state with 
respect to the goal state. We define a successful double 
check again as an infidelity below e = 10“^ and plot¬ 
ted it again as a function of the detuning A (see inset 
of Fig. [3]). This allows us to see if the approximated 
model is still good enough for optimization. Remember 
that the double check can only be lower than or equal to 
the results from RWA, since a failed optimization within 
the RWA is very unlikely to produce a successful opti¬ 
mal pulse without the RWA. Accordingly a drop in the 
double check success rate indicates the break down of 
the RWA. The behavior is similar: For small detunings, 
that is less than 10“^ we have a 100% quota of posi¬ 
tive double checks; the bigger the detuning the smaller 
the number of positive double checks. Moreover we see 
that for bigger dimensions the success rate decreases, so 
for bigger dimensions the RWA needs smaller detunings 
to still be valid. In Fig. 0] we show all optimized states 



10 ''^ 
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D 

FIG. 4: (Golor online) Reached infidelity 7 as a function 
of the distance D = \\ I’ll’a) — IV’(O)) II for different detun¬ 
ings A: 10“^ (red), 10“^ (green), 10“® (blue), 10“'^ (phik), 
10~® (cyan), 10~® (yellow), 10~^ (black). Different pictures 
correspond to different dimensions N: N = 3 (Top left), 
A = 4 (Bottom left), N = 5 (Top right), N = 6 (Bot¬ 
tom right). 

and different systems for different detunings (indicated 
by different colors, different plots correspond to different 
dimensions). On the x-axis we plot the distance between 
the initial state and the final state £* = || I'^g) ~ l'0(O)) II 
in Hilbert space, on the y-axis the infidelity / is shown. 
We see that the reached infidelity for small detunings is 
far below the threshold of e = 10“^ used for Fig. E Ad¬ 
ditionally we see the rapid drop in infidelity for smaller 
detunings and a bunching of data points around 10“^° 
or below for small detunings. In conclusion, we see that 
within the regime of a QOC task with very high fidelity 
in this important subset of possible problems can be per¬ 
formed without difficulties following the procedure laid 
out in this work. 
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V. EXAMPLE APPLICATION: ENTANGLING 
RYDBERG ATOMS 

To demonstrate the presented approach in a real phys¬ 
ical scenario we study the example of two trapped Rubid¬ 
ium atoms (23 - [^ . For each atom we consider the qubit 
states |0) = |55 'i/ 2,-F = 1), |1) = 1551/2,F’ = 2), and a 
Ryd berg: state \r) = \97D^/2) as in the experiments of 
Ref. The atoms are trapped by spatially separated 
far-off-resonance optical traps and thus can be addressed 
individually by lasers. Both qubit states can be coupled 
to the Rydberg state via two-photon transitions, e.g. via 
|5 Pi/ 2) and I 5 P 3 / 2 ). These intermediate states can be 
excluded from the model W appropriate laser detunings 
via adiabatic elimination |3l| . The atoms interact only if 
both atoms are in the Rydberg state and the interaction 
is Hint = U\rr){rr\ with U = 2'Kh ■ 20MHz, while the 
effective Rabi frequencies that couple the qubit states to 
the Rydberg state can be chosen to be a few MHz [^ . 
The system was used to experimentally implement proto¬ 
cols for CNOT-gates [H, [^, [ 2 ^ . A similar setup has 

been used to propose a protocol for multi-particle entan¬ 
glement [ 2 ^ of the type (|0...0) -I- |l...I))/-\/fV. Here we 
use the methods developed above to provide a solution 
how to transfer |00) to the Bell state (|00) -I- |ll))/-\/2 
with a single (polychromatic) pulse, so in contrast to 
most schemes we do not make use of pulsed lasers used 
sequentially. 

If we perform the RWA and individually address the 
atoms, where we drive near-resonantly the transitions to 
the four states with a single Rydberg excitation as well as 
to the doubly excited Rydberg state (we double check the 
Rydberg blockade assumption numerically), the Hamil¬ 
tonian reads 

H = ^(|00)(0r| + |lr)(10|) + ^(|00)(r0| 

+ |01)(rl|) + ^(|0r)(01| + |lr)(ll|) 

+ ^(k0)(10| -k |rl)(ll|) -k ^|0r)(rr| 

+ + F.c. 

-k 5i|0r)(0r| -k (5i -k 53)|01)(01| -k 52|r0)(r0| 

+ (^1 + <55 + 17)|rr)(rr| -k (5i -k <52 -k (53)|rl)(rl| 

-k (^2 -k (54)|I0)(I0| -k (^1 -k ^2 + <^ 4 )11'^)(1^1 

+ (5l+<52 + ^3 + ^4)|11)(11|. 

(26) 

This Hamiltonian leads to the graph that is shown in 
Fig. [5] (left). The single frequencies (same color and 
same line style in Fig. [S]) drive multiple transitions: e.g. 
the first laser (Hi, blue solid line) drives the transition 
from | 00 ) to | 0 r) as well as the transition from | 10 ) to 
jlr). As a consequence the system is not fully control¬ 
lable. However, if we switch only some of the lasers, 
namely Hi, H 4 , H 5 , and Hg, we get two controllable sub¬ 
graphs (Fig. [SJ right). One of them contains our initial 


state |00) as well as our target state (|00) -k |II))/-\/2. If 
we solve the state-to-state transfer problem using experi¬ 
mental feasible values [ 2 ^ for the Rabi frequencies as well 
as the assumption of perfect blockade with the Hamil¬ 
tonian dMl), we can achieve perfect state transfer with 
Hi = 27r/i-IMHz, H 4 = 27rF-lMHz, Hg = 27r/i-3.2MHz, 
and Hg = 2'Kh ■ 1.3 MHz, final time T = 314 ns, dg = —U 
and all other Rabi frequencies vanishing. If we double 
check the assumption of perfect blockade and go beyond 
Eq. (j26p . the infidelity increases from e = 0 to e = 0.002. 
Note that this value can be improved if single addressing 
is possible for a smaller distance between the atoms, re¬ 
sulting in a higher blockade interaction. The operation 
error is thus just a technical limitations and not intrinsic 
to the method. 



FIG. 5: (Color online) Left panel: The lasers couple the levels 
according to Eq. (I26II . Some lasers couple multiple transitions 
and the system is not fully controllable. Right panel: Using 
only four of the possible eight lasers, we are left with two 
controllable subsystems, one of them containing the initial 
state 100) as well as the state |11). On this subsystem we 
can apply the methods developed in this article to control 
the desired state-to-state transfer to the Bell state (|00) -k 

|11))/V2. 


VI. DISCUSSION AND CONCLUSIONS 

As we have seen in the previous sections QOC works 
very well within the RWA and clearly shows that the 
number of parameters for optimization has to scale lin¬ 
early with the dimension N of the Hilbert space describ¬ 
ing the system [I^ . The intuitive argument for that 

is: Assume we add an additional frequency Cj to those 
used to perform an already optimal protocol found previ¬ 
ously. If Cj is not in approximate resonance to any tran¬ 
sition frequency it can be neglected and will not affect 
the description in the RWA. If it is in approximate reso¬ 
nance to any transition frequency, we just keep the new 
frequency Cj and drop the old one which was on reso¬ 
nance with the same transition. In any case, there is no 
need to increase the number of frequencies beyond N —1. 
This supports the findings of a recent work 0 where, 
by means of an information theoretical analysis, it was 
shown that for an effective optimization the bandwidth 
of the control field (in this case given by the number of 
controls F) has to scale at least linearly with the dimen¬ 
sion of the Hilbert space associated with the optimized 
system. 
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Another connection to a recently developed and highly 
efficient optimization algorithm, namely the Chopped 
Random Basis Algorithm (CRAB) Q, can be made: In 
CRAB one expands the control pulse u{t) into a trun¬ 
cated basis, often using trigonometric polynomials with 
great effectiveness. In a typically QOC problem, one has 

77crab = Hd + u{t)Hc, (27) 

where 

Nc 

u{t) = Y. An sin {uJnt) + Bn COS {uJnt) ■ (28) 

n—1 

Analogously, one can rewrite the Hamiltonian ([T]) in the 
form 

F 

H = Hd + ^ (Re A/ cos {ojft) + Im A/ sin (w/t)) He 
/=i 

(29) 

which shows a clear one to one relationship between the 
two methods. This once again backs the observations 
made with CRAB that the number of frequencies neces¬ 
sary for good optimization results is not exceedin gly high , 
in particular it does not grow super-polynomially |l9ll20| |. 

In conclusion, we investigated the performance of QOC 
within the generalized RWA applied to a A^-dimensional 
quantum system. By introducing proper unitary trans¬ 
formations, we identified an important subset of QOC 
problems that can be described by a time-independent 
formalism, namely systems that can be described by 
acyclic, connected graphs. We solved the state-to-state 
transfer problem for the paradigmatic two-level system, 
and the dynamics of the A^-dimensional system whose 
graph is a star. We demonstrated numerically that a 
system representable by a connected, acyclic graph can 
be controlled to perform arbitrary state-to-state transfers 
and we showed that this approach allows to develop an 


optimal protocol to entangle Rydberg atoms with con¬ 
stant laser pulses, that is without the need of schemes 
for pulse shaping. Let us mention that the subset of con¬ 
nected, acyclic graphs is of natural high interest for QOC 
since they represent the class of Hamiltonians that con¬ 
nect the complete A^-dimensional system dynamics with 
the fewest possible controls, namely only A^ — 1. Re¬ 
ducing the number of controls further either leaves the 
graph unconnected or introduces a near-degeneracy, both 
impeding optimal control. 

To give an outlook, we stress that the class of QOC 
problems identified here, despite being quite general, 
does not include all scenarios where the presented 
approach can be applied successfully. In particular, 
even if the graph is not acyclic, there are cases in which 
the system can still be recast as a time-independent 
one. This is the case if the sum over all detunings 
along this cycle is zero: graphically this means that 
the phase accumulated is the same no matter which 
path of the cycle goes Q. In Ref. [I| a step further 
has been developed to encompass more of the complete 
system dynamics in a time-independent description by 
incorporating part of the counter-rotating terms, offering 
possibilities to enlarge the amount of applications even 
further. 
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